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We show that Solar System tests can place very strong constraints on K-mouflage models of gravity, which 
are coupled scalar field models with nontrivial kinetic terms that screen the fifth force in regions of large grav¬ 
itational acceleration. In particular, the bounds on the anomalous perihelion of the Moon imposes stringent 
restrictions on the K-mouflage Lagrangian density, which can be met when the contributions of higher-order op¬ 
erators in the static regime are sufficiently small. The bound on the rate of change of the gravitational strength 
in the Solar System constrains the coupling strength /3 to be smaller than 0.1. These two bounds impose tighter 
constraints than the results from the Cassini satellite and Big Bang Nucleosynthesis. Despite the Solar System 
restrictions, we show that it is possible to construct viable models with interesting cosmological predictions. 

In particular, relative to A-CDM, such models predict percent-level deviations for the clustering of matter and 
the number density of dark matter haloes. This makes these models predictive and testable by forthcoming 
observational missions. 

PACS numbers: 98.80.-k, 04.50.Kd 


I. INTRODUCTION 

K-mouflage is one of the three types of screening mecha¬ 
nisms 0] , together with chameleon and Vainshtein, for scalar 
modifications of gravity with a conformal coupling to matter 
H- Roughly speaking, these three mechanisms can be distin¬ 
guished by what triggers their implementation: K-mouflage 
operates in regions where the gravitational acceleration is 
large enough; chameleons @] are at play when the New¬ 
tonian potential is large; and Vainshtein Q5J] is active when 
the spatial curvature becomes significant. For the cases of 
K-mouflage and Vainshtein, Newtonian gravity can be pre¬ 
served deep inside the so-called K-mouflage or Vainshtein 
radii HH, which are defined as the distances from the cen¬ 
ter of a spherical source within which the nonlinearities of 
the scalar field Lagrangian become significant and, there¬ 
fore, the screening effects come into play. These two mecha¬ 
nisms share in common the fact that the screening efficiency 
(roughly determined by the size of these radii) depends solely 
on the properties of the gravitational source, as opposed to 
chameleons, for which there is also a dependence on the den¬ 
sity of the environment where the sources live. 

The existence of screening mechanisms, however, gives 
only the models a chance to pass Solar System tests of gravity 
0 ]. In other words, it is still necessary to investigate further 
the conditions that make the screening mechanisms efficient 
enough to cope with the current observational bounds. For in¬ 
stance, even well within the K-mouflage radius, the total force 
may still exhibit a non-Newtonian component (i.e. a radial de- 
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pendence that differs from 1/r 2 ), which albeit small, could 
still be large enough to induce detectable perturbations in the 
orbits of planets and moons. Currently, Lunar Laser Ranging 
experiments iHH, constrain the anomalous perihelion of the 
Moon at the 10 -11 level, which can then be used to constrain 
the K-mouflage Lagrangian density term, K(x) (cf- Eq. <[3])). 
The function K (x) is a nonlinear function (to ensure screen¬ 
ing) of the kinetic energy term \ = — [d<p) 2 /2M. i , where p 
is the K-mouflage scalar field and AT 4 the dark energy scale. 
In this paper, we shall use the anomalous perihelion of the 
Moon to determine the constraints on the shape of K (x) for 
static configurations of p. The static regime is that relevant 
for very small scales like the Solar System lfl3ll , which in our 
metric sign convention (see next section) implies \ < 0. In 
this regime, we shall also make use of the bounds imposed 
by the Cassini satellite 0 on the magnitude of fifth forces 
in the Solar System and check that we satisfy the constraints 
provided by laboratory experiments. 

On cosmological scales, the dynamics of p become impor¬ 
tant, x > 0, which broadens the range of tests that can be used 
to constrain K-mouflage models. On the one hand, in the Jor¬ 
dan frame, where matter couples to gravity minimally, New¬ 
ton’s constant becomes locally time dependent due to the cos¬ 
mological evolution of the scalar field 1 1151 1. We shall show that 
this can be used to place constraints on the coupling strength 
/3 (cf. Eq. £7}) by using the results from the same Lunar Laser 
Ranging experiment mentioned above, which currently con¬ 
strains the rate of change of the gravitational strength at the 
10 -12 yr _1 level 01 . On the other hand, the formation of 
large-scale structure is also affected by the K-mouflage field 
IU7l - [l9ll . Although the precision of cosmological data is not as 
good as that from local tests, one can still impose some con¬ 
straints by requiring that the cosmological behavior should not 
differ too much from standard A-CDM. 

In this paper, we first identify the requirements on the static 
(X < 0) and dynamical (g > 0) branches for K-mouflage 
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models to comply with the current data, and then we attempt 
to design K(\) functions that interpolate between these two 
branches in observationally and theoretically viable manners. 
A main result of this paper is that, although the small-scale 
constraints do limit significantly the functional forms allowed 
for K (x), it is nevertheless possible to construct functions that 
exhibit percent-level modifications on the growth of structure 
relative to A-CDM. This is mainly because on larger scales 
the K-mouflage screening effect becomes less efficient, which 
can have an impact on the nonlinear matter power spectrum 
and halo mass functions, as we show using the results from 
semi-analytical models of structure formation. The size of the 
differences to A-CDM are within the_ballpark of future ob¬ 
servational missions such as Euclid 0or LSST |2l. which 
makes these models predictive and therefore testable. 

The paper is organized as follows. In Sec. II, we introduce 
K-mouflage models and we summarize their main properties. 
In Sec. Ill, we focus on the small-scale regime (x < 0), which 
applies to Solar System scales, and where we consider the 
constraints from the anomalous perihelion of the Moon and 
the result from Cassini. We briefly recall the main aspects of 
the cosmological evolution and growth of large-scale structure 
in K-mouflage models in Sec. IV. In Sec. V, we analyze the 
connection between the time dependence of the gravitational 
strength and the cosmological evolution of p. In Sec. VI we 
summarize all the constraints and we build models that satisfy 
them. Sec. VII is devoted to the analysis of the cosmological 
dynamics of both the background and the matter density per¬ 
turbations in the models constructed in Sec. VI. Finally, we 
conclude in Sec. VIII. We discuss superluminality and causal¬ 
ity issues in the appendix. 


II. THE K-MOUFLAGE MODEL 


Following previous works we consider scalar 

field models where the action in the Einstein frame has the 
form 


S' = 


d 4 a 




Ml 1 


R + Rip(p) 


+ /d 4 ,v^£ m (^,^), (f) 


We use the signature (—, +, +, +) for the metric. Here, Ad 4 
is an energy scale of the order of the current energy density 
of the Universe (i.e., set by the cosmological constant), to re¬ 
cover the late-time accelerated expansion of the Universe. 

It is convenient to introduce the dimensionless kinetic en¬ 
ergy X by, 

x = ^p = a “ va ^ (4) 

Then, as described in |[T3], the canonical behavior (i.e., K ~ 
X oc —(dp) 2 / 2), with a cosmological constant pa = Ad 4 , is 
recovered at late time in the weak-x limit if we have: 

X -> 0 : K(x) - -1 + X + •••; (5) 

where the dots stand for higher-order terms. We shall impose 
this limit on all the models that we analyze. The Klein-Gordon 
equation that governs the dynamics of the scalar field p is 
obtained from the variation of the action <[T]) with respect to p 

1 H In 4 

-^=d, d»p k’] - PE = 0 , ( 6 ) 

which differs from the usual Klein-Gordon equation by hav¬ 
ing an extra term due to the coupling of the scalar field to 
matter, and where pe = —g l>v T IUJ is the Einstein-frame mat¬ 
ter density. A prime denotes partial differentiation with re¬ 
spect to x> e -g- K' = d K/dx. F° r simplicity, we assume that 
/3 = Mpidln A/dp is a constant, which implies 

A(p) = e^ v/Mpl . (7) 

The normalization of the first two terms in Eq.© only de¬ 
fines the normalization of the constant Ad 4 and of the field p 
(except for its sign), and therefore it does not imply any loss of 
generality. Similarly, the sign of (3 in Eq.© can be absorbed 
in the sign of p, and therefore we can choose ,7 > 0 without 
any loss of generality. 

For observationally interesting cases, we have \A— 1| < 0.1 
(see for instance lfl7l4l9ll or the discussion in Sec.[V]). There¬ 
fore, higher-order terms in In A(p) would only make small 
quantitative changes with a negligible impact on our conclu¬ 
sions. Next, we shall first analyze the solutions to the Klein- 
Gordon equation for spherically symmetric static configura¬ 
tions. Then we investigate the solutions in a cosmological 
setting, where the dynamics of p become important. 


where g is the determinant of the Einstein-frame metric tensor 
and Urn ' 1 are various matter fields defined in the Jordan 
frame. This also defines the Einstein-frame Newton’s constant 
as 87 tQ = Mp 2 . The K-mouflage scalar field p is explicitly 
coupled to matter through the Jordan-frame metric g^ v , which 
is given by the conformal rescaling 

5 ^ 1 / = -A (p) g^is, ( 2 ) 

and g is its determinant (here and throughout, quantities with 
a tilde are defined in the Jordan frame). In this paper, we 
consider models with a non-standard kinetic term 

C v (p) = M 4 K with X =-^d^pdu.p. (3) 


III. LOCAL DYNAMICS 

In this section we study the model predictions for Solar Sys¬ 
tem scales, where we work in the Einstein frame. We assume 
the fields are static, which is a valid approximation since we 
consider time scales that are much shorter than the cosmolog¬ 
ical ones. In this case, A(p), associated with the slow run¬ 
ning of the cosmological background, is approximately con¬ 
stant and the Einstein and Jordan frames are equivalent. 1 To 


1 In Sec.[V] we shall relax the static approximation when we study the slow 
time variation of Newton’s gravitational strength in the Jordan frame. 
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simplify notations we drop the subscript “E” in this section 
and we denote c as the speed of light. One can split the scalar 
field as 


¥>iocai(r, t) = p(t) + p( r), (8) 

where p is the uniform value associated with the cosmological 
background and p( r) is the perturbed component on which we 
focus in this section. The cosmological background plays no 
role in this section, apart from setting the value of the coupling 
factor (3 ~ d\nA/dp. 


A. Static case 

For a source with density p, the static Klein-Gordon equa¬ 
tion becomes (see III for details) 

= (9) 

with \ = — c 2 * (V r </?) 2 /(2.M 4 ), from which one can obtain a 
first-order algebraic equation for \7 r p 

v rV r = ^(v r % + v r xw), go) 

Cr 

where T'n is the Newtonian potential, given by the usual Pois¬ 
son equation 


V^ n = 4tt Gp, ( 11 ) 

and u is a divergence-free potential vector (which must be 
determined along with p). 

The right-hand side of the Poisson equation also involves 
the fluctuations of the scalar field energy density 5p v , but as 
shown in Ref. Q, this would only introduce negligible ef¬ 
fects compared to the fluctuations of the matter density. 

The spatial fluctuations of the coupling function A(p) can 
also be neglected in most cases (see fl3lll8ll and Eq. ( l54l > be¬ 
low), except in the Euler equation or the geodesic equation 
below, which involve the gradient of A{p) that gives rise 
to the fifth force associated with the scalar field gradient. In 
a similar fashion, we work in the weak gravitational field and 
non-relativistic limit, so that the metric fluctuations only 
appear in the Euler equation or the geodesic equation (fl9] > be¬ 
low through the gradient V r \P]sr- For the spherical configura¬ 
tions we consider here, Co = 0, which allows one to analyze 
the dynamics of the system and the scalar force due to p in a 
more straightforward manner. 

The fifth force generated by the K-mouflage field is given 
by El 


/?c 2 2/3 2 

p ^-fe v ^ = -Tr v '* N - (12) 

The K-mouflage screening mechanism relies on the fact that, 
in the nonlinear regime, i.e. deep inside the K-mouflage ra¬ 
dius, the factor K' becomes large, which suppresses the fifth 
force relative to the Newtonian one, Fn = -V r \F;N (note that 


IF^I ~ |Fn/A"'|). The implementation of the screening can 
be illustrated in a few steps (see also (H). For a spherical 
matter distribution p(r) with mass profile 

M(r) = f dr' 47 rr ' 2 p(r'), (13) 

J o 

the Klein-Gordon equation i flOli can be written as 


dp , _ f3M(r) 
dr c 2 Mp\Airr 2 ’ 


with x = 


2 M 4 



(14) 


We define the “K-mouflage screening radius” Rx by El El 


Rk 


{ /3M \ 1/2 

\4ircM Pl M 2 J ’ 


(15) 


where M = M(R) is the total mass of the object of radius 
R. Then, by introducing the rescaled dimensionless variables 

x = t/Rk, m(x) = M(r)/M and 4>{x) = <p(r)/ipK, with 

Pk = M 2 Rk/c, (16) 

the integrated Klein-Gordon equation (IT 4 l) becomes 




with x 


1 ( 

2 \dxJ ’ 


(17) 


which can also be written as 

P^K'bc) = 08 ) 

X 2 

A unique solution of the Klein-Gordon equation is always 
guaranteed when y/—2xK'(x ) i s a monotonic decreasing 
function over x < 0 , which grows up to -foo as x ~oo. 
This is assumed in the following 2 and it implies K'+2xK" > 
0. A test particle outside the matter distribution evolves ac¬ 
cording to the non-relativistic equation of motion 

d 2 r Be 2, 

d? = (19) 

which becomes the same as in standard gravity, provided one 
interprets the equation with a total potential that is the sum of 
the fifth force one d T 

Be 2 

^ (20) 

M pi 

with the Newtonian potential 'Fn- For a spherical body we 
can consider radial trajectories and the scalar field gradient is 
given by Eq. O 

dp = p K m(x) 
dr Rk x 2 K r 


2 See IfTl for a more detailed study of this case, and the relaxation of the 
scalar field to its equilibrium state, as well as a discussion of badly be¬ 

haved cases where there are no well-defined static profiles, i.e. when 
x/—2xK' (x) is n ot a monotonic function that decreases from +00 to zero, 

as x varies from —00 to zero. 
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Outside the spherical body we have m(x) = 1 and we obtain 
the equation of motion 


dV__£M / 2/3 2 \ 

df 2 r 2 \ K'(x(r)) J 


( 22 ) 


This equation shows that the standard gravitational law ac¬ 
quires extra terms proportional to /3 2 /K'. This can be cap¬ 
tured by defining an effective gravitational strength 


G eS (r) 



W 2 \ 
K '(x(r))J 


G 


(23) 


which depends on the distance from the central object. As the 
test particle gets deeper inside the K-mouflage radius, then 
m/x 2 3> 1 and K'(x( r )) becomes larger. If K'(x( r )) be¬ 
comes large enough, then the correction term in Eq. (l23l) be¬ 
comes negligible and Newton’s law is retrieved. Next, we in¬ 
vestigate the conditions for this correction to be small enough 
to evade the constraints on the anomalous perihelion of the 
Moon set by Lunar Laser Ranging experiments lf22j] . 


B. Perihelion constraints 


1. Constraint on the kinetic function 


The explicit dependence on r of the correction to Newton’s 
law implies that orbits are perturbed and in particular that 
there is an anomalous perihelion for objects like the Moon 
orbiting around the Earth. The perihelion angular advance is 
given by 
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(24) 


then, its derivative with respect to \ gives the spatial variation 
of x as 


da; x c 2 

d\ 4x c 2 ’ 


(30) 


where we have defined the jmeed of scalar spherical waves 
around the massive object 111311 


c? = 


K' + 2 xK’ 
~K’ 


2 ^ 2 
■C > C . 


(31) 


This yields our final result for the anomalous perihelion ad¬ 
vance 


69 = — 8-7T 


/ 3 2 XK" c 2 
K’ K’ cV 


(32) 


as expressed solely in terms of x and of the Lagrangian func¬ 
tion K and its derivatives. This expression contains the factor 
/3 2 /K', which, as we have seen in Eq. (l23l > controls the ampli¬ 
tude of the fifth force. The Cassini satellite results constrain 
this amplitude in the Solar System in such a way that < 
10" 5 GJ]. On the other hand, the Lunar Ranging experiment 
implies for the Earth-Moon system that \60\ < 2 x 10 -11 . 
Thus, the Lunar Ranging constraint is much stronger than the 
Cassini bound and will prove to be the main source of con¬ 
straints on the form of K (x) on the static branch X < 0 
(cf. Sec. lV’ll below). The only sensible way of complying with 
this bound is to reduce xK"/K’ in the static case, which also 
gives c s ~ c. This can be achieved by suppressing the contri¬ 
bution of the nonlinear terms in K(x) when x < 0. 


2. Constraint on the running of the coupling /3 


where e is the ratio between the fifth-force and Newtonian 
potentials 


(iT (3c 2 tp (3c 2 rp 

4/n .Mpi'I'n MpiQM 


(25) 


and where we have used T'n = — GM/r. This gives explicitly 

£ f3c 2 ip 


r MpiGM ’ 

and, using Eq. (fl4l) . one gets the variation 

2/3 2 


_d 

dr 


(?) 


K’r 2 


This implies that the anomalous perihelion is given by 


59 = 7rr— 
dr 


— 2/9 2 

K' 


- 9 ^ 

2 ’ K’ 2 K Ax' 


(26) 


(27) 


(28) 


where, we recall, x = t/Rk■ Using the fact that, outside the 
spherical source the Klein-Gordon equation (fill) yields 


V^xK'(x) = ^ 



(29) 


Throughout this paper, we focus on the case of a constant 
coupling strength j3, which corresponds to the exponential 
coupling function £7]). However, we may also consider more 
general coupling functions where /3 would now depend on 
time and space through the variations of the scalar field ip. 
This would not change our results for the kinetic function 
K(x) and the typical amplitude of the coupling /3, because 
\ip/Mp\\ (which goes to zero at high redshift) does not grow 
beyond 0.5 for observationally interesting models, as seen in 
Eqs. (l48l) and (l54l > below. Hence the variations of /3 are small 
in realistic models. 

On the other hand, we can investigate whether the very 
small bound on the anomalous perihelion, \50\ < 2 x 10 -11 , 
provides interesting constraints on the possible amount of run¬ 
ning of /3. For a general coupling function A(p) the fifth-force 
potential reads as 6 T = c 2 hi(.4/.4) and the perihelion ad¬ 
vance dM gives 
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r 2 c 2 /3 dp 
GM Mp\ dr 


(33) 


where we used again 4r N = —GM/r. If /3 depends on space, 
the Klein-Gordon equation ((9]i cannot be exactly integrated as 
in Eq. dTfll . However, the fluctuations of i3 can be neglected 
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at first order, and hence we can still approximate dip/dr by 
Eq.(fl4l). This yields 
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(3 d/3 /3 2 d K' 

+ (34) 


The second term is the one that was already obtained in 
Ea. (l28l) for constant /3. Focusing on the first term and using 
again Eq.i fl4l i. we obtain the contribution to the anomalous 
perihelion due to the spatial variation of /? as 


68 = 8n 


*n P 2 P' 

c 2 K' 2 ' 


(35) 


Here we defined 


P 


M P \ 


din A 

dip 


P' 


Mpi 


dp 

dip 


— Mp j 


d 2 In A 
dp 2 


(36) 


the dimensionless derivative of the coupling with respect to 
ip/Mp\. We shall see in Sec. IVICl that. because of the small 
prefactor /c 2 , the bound on the perihelion advance only 
gives a very loose bound on the derivative /3', which does 
not provide useful information on the coupling function A{ip). 
This is consistent with the fact that fluctuations of the scalar 
field ip and of A can be neglected in most places, except as the 
source of the fifth force that explicitly involves the gradient of 
A. 


C. Laboratory tests 

Measurements of the orbits of planets in the Solar Sys¬ 
tem constrain the deviations from General Relativity, through 
Eq.(l23l) or Ea. (l32l) . In addition, laboratory experiments, such 
as the ones using atom interferometry and measuring the ac¬ 
celeration induced by a test mass of a few grams over dis¬ 
tances of a few centimeters, also constrain the amplitude of 
the fifth force in Eq. (l23l) to a 10 -4 accuracy [23i, 24]. Thus 
laboratory experiments place constraints on 2 /3 2 / K' in the 
static case, but in an even more non-linear regime than the 
Cassini spacecraft or the Lunar Ranging experiment. Indeed, 
laboratories on the surface of the Earth are further inside the 
K-mouflage radius of the Earth than the Moon. 


IV. COSMOLOGY 

Before analyzing the constraints on K-mouflage models ob¬ 
tained from regimes where the dynamics of ip are important 
(X > 0), we first briefly recap the equations relevant for the 
cosmological evolution of the background and of linear per¬ 
turbations. For further details, we refer the reader to JH for a 
study of the background expansion history, to [; lj| for a study 
of large-scale structure formation and to JH for a study of 
the model predictions for the Cosmic Microwave Background 
(CMB) temperature and lensing potential power spectra. 


A. Background 


Considering only spatially-flat universes, the Einstein equa¬ 
tions lead to the usual Friedmann equations [17]], 

3MfH 2 = Pe + Pp, (37) 

—2Mp\H = pe + Pip + pip (38) 

where pE, Pp and p v , are, respectively, the matter and scalar 
field energy densities and pressure in the Einstein frame: 

p v = -M 4 K + p 2 K’, pip = M 4 K. (39) 


The overbar denotes uniform background quantities, and the 
dimensionless field \ (Eq. (|4}) simplifies as 


X = 



(40) 


It is convenient to introduce the rescaled matter density p. 


p = A 4 p E , 


(41) 


which satisfies the standard conservation equation. The Klein- 
Gordon equation <[6]) is now given by 

5 («■#*•>-$*■■• (42) 

We can also define an effective scalar field energy density 

pf = Pip + [A{ip) - 1 ]p, (43) 

which satisfies the standard conservation equation (the pres¬ 
sure pip is not modified Il7jl) 

pf = -3 H(pf + Pip). (44) 

Under these definitions, the Friedmann equations (l37l i- (l38ll 
become 


3M^H 2 = p + pf.. (45) 

-2 M 2 x H = p + pf+p v . (46) 

We define also the time-dependent cosmological parameters 


— 


p + pf ’ 


= 


Pf 


p+p: 


;eff > 


w eS = 

Hip 


(47) 


At early times, we have </7 —y 0 and A[ip) is normalized by 
aw = i m. For observationally interesting models, we 
have A ~ 1 (|A — 1| < 0.1) until today fl7l[l9ll (see also the 
discussion on the constraints from Big Bang Nucleosynthesis 
(BBN) below). From Eq. (l42l >. one can then write 


Ppt_ Pv _ Pf 

Mp\K' ’ Mpi ~ K' ‘ 


( 48 ) 
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B. Linear Perturbations 


On large scales, the evolution modes, D +(n) , of small lin¬ 
ear density fluctuations satisfy the equation 11811 


d 2 P 

dr] 2 


1 - 3 wfflf 


f-2 | —— —f2 m (l + ei)D — 0, 


(49) 

where 77 = In a is the number of e-foldings. The functions 
and £2 are time dependent only and they are given by 


and 


i 11 2 ^ 2 


P 2 

" ' K' 


K' 


d In A f} d tp /? 2 

dr] Mpi dr] K' 


(50) 


(51) 


In Eq. (f50t the sign of e 1 cannot be determined a priori be¬ 
cause the terms {A — 1) ~ /3<p/Mp\ and 2Af3 2 jK' are typ¬ 
ically of opposite signs and of the same order. In Eq. ( |49| . 
ei mimics the effects of a modified Newton’s gravitational 
strength on the linear growth of structure, while £2 appears 
as a friction term. Both these terms are of order (3 2 /K', just 
like the case of (p in Eq. ([48}. 

The reader might note that whereas in the cosmological lin¬ 
ear equation ( 149} Newton’s gravitational strength appears to 
be amplified by a factor (1 + ei), in the static case (123} we 
found a factor (1 + 2(3 2 /K') and in Eq. (f52l) below we have 
a slow drift given by A 2 . These differences come from the 
fact that they apply to different regimes, which also involve 
different sets of approximations. This explicitly shows that 
modifications of gravity, as defined from the Lagrangian of 
the theory, can have subtle effects that are not captured by a 
unique rescaling of Newton’s gravitational strength. 

More precisely, the linear evolution equation (149} applies 
to the cosmological density field in the Einstein frame, where 
the energy-momentum tensor is not conserved because of the 
time-dependent conformal mapping ([2} to the Jordan frame 
that defines the matter Lagrangian. This slow drift between 
the Einstein and Jordan frames gives rise to the factors A in 
Eqs.(|5Q} and (BT} . In contrast, the static equation (123} ap¬ 
plies to small-scale systems over time scales that are short as 
compared with the Hubble time, so that the cosmological vari¬ 
ation of A and the expansion of the Universe can be neglected. 
Then, non-conservation terms of the form pd In A/dt can be 
neglected and there is no difference between the Einstein- and 
Jordan-frame density fields [thus, with A = 1 Eq.( |5 11 ) gives 
£2 = 0 and Eq. ([50} gives back £1 = 1+2 j3 2 /K’, as in Eq. ([23} . 
except that K 1 is the cosmological background value of the 
kinetic function whereas in Eq.(l23} it is the local perturbed 
value associated with the compact object]. Finally, Eq.(l52t 
gives the cosmological drift of Newton’s coupling in the Jor¬ 
dan frame, as opposed to the Einstein-frame evolution equa¬ 
tion ([49} (where by definition Newton’s coupling is indeed 
constant) and to the static equation ( [23} (where cosmological 
drifts are neglected as compared with the system time scale. 


such as the orbital period of the planet). Thus, these differ¬ 
ences come from the fact that we consider different time and 
length scales, and different frames. 


V. TIME VARIATION OF Q 


In Sec. m we have neglected the effects of the cosmolog¬ 
ical evolution on the predictions for the Solar System, since 
we considered time scales that are much shorter than the Hub¬ 
ble scale. However, in K-mouflage models, the conformal 
mapping of Eq. ([2} implies that in the Jordan frame, where 
matter couples minimally, Newton’s gravitational strength, Q, 
should be time varying. This can be understood as follows. 
Let us consider the conformal transformation from the Ein¬ 
stein to the Jordan frame. Since, the Ricci scalar and ■sj -g 
transform, respectively, as R = A 2 R + ..., (where the dots 
stand for additional terms associated with derivatives of A), 
and sj—g = A~ A y/—g, then one has Q = A 2 Q = A 2 /8nMpp 
That is, Q becomes time-varying due to the background evo¬ 
lution of A{(p). This time variation of Newton’s gravitational 
strength in the Jordan frame can be constrained in two ways: 
(a) through the comparison between the local value of Q and 
that at the time of BBN, and (b) through the impact on the 
trajectories of planets and moons in the Solar System. 

Starting with the BBN case, a value of Newton’s constant 
which would be different during BBN compared to the one 
inferred from local measurements would be tantamount to a 
change in the Hubble rate and therefore would lead to a dis¬ 
crepancy in the formation of the elements. Such a change 
cannot exceed about ten percent [25J 26j|. In the Jordan frame 
we have 


which implies the bound 

■77 “ l+BBN — +local| ^ 0.05. 
lVlp\ 


(52) 


(53) 


At the time of BBN, we have +bbn = p{z ~ 10 10 ) ~ 0. The 
local value of the scalar field is given by Eq.®. In Sec. HIH 
we neglected the contribution from the background part, (p, 
since we were interested only in the additional ?’-dependence 
in the force law. Here, however, we should take it into ac¬ 
count, and from Eq. ( 148} we have (3(p/Mp\ ~ (3 2 / K’ ~ /3 2 
(because K' ~ 1 at low redshifts). From Eq. (IT2} . we have 
that the perturbed part of the scalar field, tp(r), is of order 
/3p(r)/Mpi ~ /3 2 '&n/K'c 2 -C /3 2 , because K' 1 in the 
Solar System (see Sec. IVIAl below) and T'n/c 2 <C 1 (in the 
Solar System we have 4 /n/c 2 ~ 10~ 6 ). Altogether we have 


|+(r)|-C |+(f)| and A(^i oca i) ~ A, (54) 


and Eq. ( [53} implies 


/3 2 < 0.05. 


(55) 


A second type of constraint on the rate of change of the 
gravitational strength comes from the change with time of the 
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trajectories of planets and moons. This has been monitored 
by the Lunar Laser Ranging experiments for the Earth-Moon 
system 0. Recalling that Q ~ A 2 Q, then its rate of change 
is given by 


ding 

df 


2e 2 
1 + e 2 


H 


2 e 2 H, 


(56) 


VI. CONSTRAINING THE MODELS 


In this section, we summarize the constraints on K- 
mouflage models discussed above and also other stability con¬ 
ditions presented previously in the literature iIdL 1171 11811 We 
also build models that satisfy them. 


where H and t are the Hubble rate and the time in the Jor¬ 
dan frame (which are related to the Einstein-frame ones by 
H = H(l+e 2 )/A and df = Adt), and we have used Eq. ( BTI) . 
Hence we find a direct link between the behavior of cosmolog¬ 
ical density perturbations and the Lunar Ranging constraint. 
The current bound gives that |dlng/df| now < 1CP 12 yr _1 
lfl6tl . In particular, taking as a reference value h = 0.67 (al¬ 
though this value is not critical for the conclusions) we find 
that the constraint on g 2 reads 

M now < 7.3 x 10- 3 . (57) 

From Eq. ( |5H we can see that this gives a constraint on the 
ratio p 2 /K' today. In fact, this is a strong constraint on the 
coupling P of the K-mouflage field to matter, which is inde¬ 
pendent of the details of the kinetic function K(\)- At late 
times, x g° es to zero, and we typically have x <C 1 today, as 
well as, K' ~ 1, from Eq. (01. Therefore, we find 

£<10-\ (58) 

which is a tighter bound than the BBN constraint of Eq. (l55l >. 

In the Einstein frame, the variation of Q with time is re¬ 
placed by a variation of the masses of fundamental particles 
as = Am ^ where m 3 is the mass in the Jordan frame. 
In particular, we can see that the ratio of fundamen¬ 

tal particle masses over the Planck scale is frame-invariant. 
For bound states such as the protons, and as long as the QCD 
phase transition can be modeled in the Jordan frame where no 
coupling of gluons to the scalar field are present, the masses in 
the Einstein frame are still proportional to A. This guarantees 
that the proton-to-electron mass ratio m p /m e is independent 
of time and no constraint on K-mouflage can be drawn from 
the bounds on the variations of rn p /m e from quasar absorp¬ 
tion lines. 

From Eqs. d50| and (ED, ones notes that ei and e 2 are of 
the same order. Consequently, the bound that Lunar Rang¬ 
ing tests place on e 2 , translates also into similar bounds to 
6\, which highlights an interesting connection between Solar 
System constraints and the growth of structure on cosmolog¬ 
ical scales. In particular, from the above bounds one expects 
deviations from standard A-CDM to be of the order of a few 
percent 00. On the other hand, the Lunar Ranging test 
does not constrain the past behavior of the model, when x 
deviates from zero. Therefore one must still investigate the 
cosmological impact of the functional form of K(x) when 
X > 1. This shall be done in part of the discussion of the next 
section. 


A. Combined Constraints 


The cosmological regime corresponds to \ > 0, with 
X Too at early times and x —> 0 at late times. To avoid 
ghosts we must have K' > 0 and K' T 2 xK" > 0 (l3], as 
fluctuations Sip around the cosmological background <p prop¬ 
agate with a speed c s given by 


c 2 = 


K' 


K' T 2 xK" 


2/2 
-C < C . 


(59) 


This is formally the inverse of Eo. (l3T1) . but in Eo.(l59l> x — 
X > 0 is the homogeneous time-dependent cosmological 
background, whereas in Eq.d59t < 0 is the small-scale static 
solution. The function y/x,K f (x) mus t also increase monoton- 
ically up to Too, so that we have a well-defined cosmological 
behavior up to high redshifts ii, which again implies 
K' T 2 \K” > 0. Requiring that the dark energy component 
becomes negligible with respect to the matter density at early 
times implies that I\(x) grows faster than x (e.g., as a power 
law x m with m > 1, see Ref. J0). The marginal case, where 
p v grows as 1 / 1 2 in the early matter era (like the matter den¬ 
sity) but is a small fraction of the matter density, corresponds 
to a constant K' with K' T> 1. Moreover, one can also impose 
that the growth of large-scale structures should not differ too 
much from A-CDM. A reasonable maximum deviation can 
be placed at the few-percent level today, which leads to the 
bound P 2 /K' < HP 2 , from Eqs. (|50|)-(|5I1). This requires 
p 2 < 10 2 because K' ~ 1 today. This bound also ensures 
that the BBN constraint, Eq. (l55l) . is satisfied. Therefore, the 
cosmological constraints are: 


cosmological constraints: semi-axis x > 0, 


K' >0, K' T 2 X K" > 0, (60) 


V\K'(x) -A Too for x -A Too, (61) 

K' > 1 for x > 1, /3 < 0.1. (62) 

The small-scale static regime corresponds to x < 0. To 
avoid singular behaviors and to ensure well-defined solutions 
for any matter density profiles, we must have K' > 0 and 
\J~X_A'(x) monotonic and unbounded over x < 0 Q [re¬ 
call also the discussion about Eq. (IT8l )l. This latter condition 
also implies K' + 2 xK" > 0 and ensures that the c 2 > 0, 
where c s given in Eq. ( 13 1 1 ) is the propagation speed around a 
static background. Thus, we have the very general small-scale 
conditions: 

small-scale constraints: semi-axis x < 0, 
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K' > 0, K' + 2 X K" > 0, (63) 

V~X K '(x) -> +oo for x —» -oo. (64) 

Note that instead of \f~~xK'ix) +°° at infinity it may be 
sufficient to go to infinity at a finite negative value of In ad¬ 
dition, the Solar System dynamics corresponds to x < —X*> 
as we are far in the nonlinear regime, where we denote by 
(—X*) the transition between the linear regime, where I\' ~ 1 
and the kinetic function is governed by the first two terms of 
the expansion <(5]), and the nonlinear regime where K' 1: 

small-scale highly nonlinear regime: \ < — x*, 
with — x* < 0 and K'(—x*) 3> 1. (65) 

The Cassini bound on the amplitude of fifth forces in the 
Solar System implies that the scalar field should be suffi¬ 
ciently screened locally. From Eq. ( 1231) this requires that 

/3 2 

Solar System screening: — < 10 5 for x ~ Xs.s. (66) 
A ' 

where K' must be evaluated at values Xs.s. that correspond to 
the Solar System regime, that is, at distances of order one AU 
from the Sun. This is well within the K-mouflage radius of the 
Sun itself, which means that Xs.s. < — X* an< i large values of 
K’ in this regime are consistent with the requirement K'( 0) = 
1. More precisely, using that A4 4 ~ pdeo is roughly the mean 
dark energy density today, then we obtain from Eq. (IT5l) that 
the K-mouflage radius of an object of mass M is given by 


Rk{M) ~ J — 3470 AU. (67) 

V 

Thus, for f3 ~ 0.1, the K-mouflage radius of the Sun is 
Rk(Mq) ~ 1097 AU, which is much larger than the size of 
the orbits of all Solar System major bodies (Neptune and Pluto 
being at about 30 and 40 AU). Moreover, the integrated Klein- 
Gordon equation ( l29l ) gives for the Solar System regime: 

V=X^.K. S . ~ 10 6 . ( 68 ) 

In practice, the constraint (l66t means that we require 

K' » 1 for x « -X*, (69) 

which automatically ensures that the general conditions ( l64l i 
are satisfied. Alternatively, one could also have very small 
values of the coupling parameter /3, but this would yield a 
cosmology virtually indistinguishable from A-CDM, render¬ 
ing the K-mouflage scenario less interesting. In particular, if 
we assume (3 ~ 0.1, so that deviations from the A-CDM cos¬ 
mology are not completely negligible, we require K' > 10 3 
for X < -X*- 

From the constraints on the anomalous perihelion of the 
Moon, Eq. (l32l) . we have 

Moon perihelion: --— < g x 10” 13 

F K' K' + 2 X K" ~ 

for X~Xm.e., 


where Xm.e. corresponds to the Earth-Moon system. The K- 
mouflage radius of the Earth is about 2 AU, and the distance 
between the Earth and the Moon is d ~ 2.6 x 10 -3 AU, which 
gives [Rk(M^)/ d] 2 ~ 0.6 x 10 6 . Therefore, the Sun and the 
Earth have about the same impact on the scalar field config¬ 
uration at the location of the Moon. In practice, this means 
that the value of x associated with the perihelion constraint 
is roughly the same as that associated with the Cassini bound 
©: 

Xm.e. ~ Xs.s.- (71) 

In other words, the constraints of Eqs. ( l66l > and (f70t both apply 
to the regime set by Eq. ( l68l >. 

Laboratory experiments give a constraint on the Newtonian 
force of order 10 -4 , which means from Eq.( l23| i that we have 

3 2 

laboratory: — < 10 -4 for x ~ Xlab- (72) 
A 

In this case, the screening is induced by the gravitational field 
of the Earth. Since Rk(M^) ~ 2 AU, we have at the surface 
of the Earth [Rk (M®)/f?®] 2 ~2x 10 9 , and Eq.(l29l) gives 

- 10 9 - (73) 

This means that laboratory experiments constrain the K- 
mouflage model much further into the nonlinear regime than 
the Cassini or Lunar Ranging probes, with typically | Xlab | 2^ 

I Xs.s. | • As described in Sec. IVIIAl below, as explicit mod¬ 
els that pass all constraints we shall consider simple models 
where K' converges to a large constant value K t in the non¬ 
linear regime |x| X*- Then, the constraint d72t on K' is 

less stringent but close to the Solar System one (l66t . 

Finally, the Lunar Laser Ranging constraint on the local rate 
of change of Newton’s gravitational strength, Eq. (l58l) . implies 

bound on the time dependence of Q : (3 < 0.1. (74) 

B. Possible kinetic functions 

From Eqs. (l62l) and ([74]), we note that the cosmological and 
local (Earth-Moon system) constraints on the coupling param¬ 
eter happen to be of the same order, 3 0.1. In terms of 

the kinetic function K(x), the cosmological and small-scale 
constraints apply to different branches, x > 0 and x < 0, re¬ 
spectively. Therefore, there seems to remain some freedom in 
the choice of the function K. The main requirements are that 
K' should be large in both limits x -7► ±oo, which ensures 
screening in both the early-time cosmology and the small- 
scale dynamics. However, in addition to this, we also have 
that in the small-scale regime, around Xs.s. < 0, the kinetic 
function is very strongly constrained by the perihelion bound, 
Eq. (®. 

In order for the function K (x) to satisfy the above require¬ 
ments, one cannot avoid some degree of fine-tuning. A sim¬ 
ple way to see this is to note that power-law behaviors cannot 
easily match the constraints. For instance, the Solar System 
(70) regime, Eq. (l68l >. requires a large degree of nonlinearity for 
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K(x), far away from the low-order expansion of Eq. [5] which 
would give K' ~ 1 and (—Xs.s.) ~ 10 12 from Eq. (l68l) . and 
would fail to satisfy the screening criterion, Eq. {66}. This 
suggests that (—Xs.s.) 7?> 1 is far in the nonlinear regime 
of the function K(\), where K' is also much greater than 
unity. However, this is not sufficient to satisfy the perihelion 
constraint, Eq. ( ITOl ). Thus, considering a power-law behav¬ 
ior K'{x) ~ ( — x) m > with m > 0 , we have |xiT"| ~ K'. 
The perihelion constraint, Eq. ( 1701 . is much stronger than the 
Cassini constraint, Eq. {66} , and we obtain for 3 ~ 0.1 the 
lower bound K' > 10 10 . Then, Eq. {68} would actually imply 
( — Xs.s.) < 10^ 8 . Therefore, such a power-law regime would 
need to occur very close to the origin, with K' going from 1 
to 10 10 as x goes from 0 to —10 -8 , that is, x* < 10 -8 . This 
would be an extremely finely tuned kinetic function K (x). 

To achieve this quick growth of A' 7 (x), we may consider 
functions that diverge at a finite negative value —x* < 0, such 
as K'(\ ;) ~ (x + X*)~ m with m > 0. Then, saturating the 
upper bound of Eq. {70} with the condition of Eq. {68} gives 
(x + X*) ~ 10 -10 (for m = 1) and x* ~ 10 -8 . There¬ 
fore, including singular kinetic functions does not remove the 
need for extreme fine-tuning and again requires a very quick 
departure from the low-x expansion ©. 

The way out is to suppress the second derivative K" , that 
is, to look for kinetic functions such that 

X«-X*: \xK”\ « K' . (75) 

This means that ln(AT') must grow much more slowly than 
ln(—x) for x -A ~ oo. Typically, this corresponds to mod¬ 
els where K' converges to a constant (although a logarithmic 
growth could also be possible), 

X —> —oo : K' -A if* > 1, (76) 

with the relation {68} giving K t ~ 10 6 / y/— Xs.s.- Then, the 
Cassini bound, Eq. {66} . implies for 3 ~ 0.1: 


between 10 3 and 10 6 at large negative x- This transition, 
however, does not need to be very sharp, as it can take place 
over an interval that can be as large as x* ~ 10 6 . Never¬ 
theless, it still requires introducing a parameter K t > 10 3 
for the asymptotic constant slope of the kinetic function. We 
note also that this cannot be obtained using a logarithmic 
growth of K', which, although consistent with Eq. {75}. as 
ln(x*) < ln(10 6 ), it is not sufficient to generate factors of 
order 10 3 . 


C. Running of the coupling 3 

The Newtonian potential of the Sun at the orbit of the Earth 
is T'n/c 2 — -10 -8 (the potential of the Earth at the orbit 
of the Moon gives the smaller value — 10~ n ). Then, Eg. (031) 
gives the constraint 

\3'\ < 8 x KT 5 . (80) 

For 3 = 0-1 and K * = 10 3 this yields 

\3'\ < 8000. (81) 

Therefore, in contrast with the kinetic function K{x), the 
bound on the perihelion advance does not provide useful con¬ 
straints on the shape of the coupling function A(<p) [generic 
functions of the form A(j3(p/Mp{) would have & ~ /3 2 ~ 
10 2 ]. This is due to the small prefactor \I7 n/c 2 that appears 
in Ea.(l35}. 

VII. EXPLICIT MODELS THAT PASS ALL 
CONSTRAINTS 

A. Constructing models 


if 3 - 0.1 : fC* > 10 3 and x* < 10 6 , (77) 

which also ensures that the laboratory constraint d72} is satis¬ 
fied. On the other hand, if we require (— Xs.s.) 3i 1> so that 
the transition from K'( 0) = 1 to K, does not take place in a 
very small interval, to avoid extreme fine-tuning, we have the 
upper bound 

if X* > 1 : K, < 10 6 . (78) 


A family of models which satisfy the properties obtained in 
Sec. lVIBl is given by 

n 

K' = 1 + iv*—--, n even integer, (82) 

X* + X n 

which is well defined over the full real axis (hence the choice 
of even integers for the exponent n ) and goes to at large 
|x|, with the constraints 


Thus, we obtain a finite range for the possible values of iC*. In 
particular, it happens that too large values of K' are ruled out 
if we wish to avoid too much fine-tuning. Moreover, Eq. {70} 
yields 


IXs.s.A^J < 1 Q-10 Jsr / s < 10 -4 J ( 79) 

^S.s. 

which means that we must have converged to the asymptotic 
value of K' at the 10 -4 level at least. 

Admittedly, there still remains some tuning, as we require 
that K' (x) g° es from unity at x = 0 to an asymptotic value 


10 12 

(9 = 0.1, A*>10 3 , X*<-F5- (83) 

K; 

This ensures that the Cassini bound (166} is satisfied with 
I Xs.s. | > X* from Eq. (l68k as well as the laboratory constraint 
{72} . Then, the perihelion constraint {70} yields 


X* < 



10 12 
Ki ' 


(84) 


As seen in Eq. {78} . we should have A'* < 10 6 if we wish 
to avoid tuning x* to values that are smaller than unity. This 
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FIG. 1: Upper panel: kinetic function K{\) for the arctan 
models of Eq. 1861 . with (A'*,X*) = (10 3 ,1) (I, blue and 

crosses), (A*,x») = (10 4 , 1) (II, brown and squares), (K t . x*) = 
(10 3 ,10 2 ) (III, green and circles), and the cubic model of Eq. ( I94l > 
with Ko = 1 (IV, red and continuous line). We show 1 + K, which 
is an odd function of x, with 1 + K > 0 for x > 0. Lower panel: 
derivative K'(x) f° r the same models. It is an even function of x- 


agrees with Eq. (l84l) . which in such cases is more stringent 
than the last equation in (l83l> . In particular, this gives: 

10 7 10 12 
n = 2: X* <-^ 3 / 2 > n^-oo: X* < -Jq- (85) 

In the remainder of this section we focus on the simplest 
model with n — 2, which corresponds to 


K( X ) = -1 + X + K*[X - X* arctan(x/x*)], (86) 


with the low-x expansion 

X-fO: K(x) = - l + X+^^2 

3 xi 


5 xi 


We consider the following three models. 


(I): 

A'» 

= 10 3 , 

X* 

= 1, 

(II): 

AT* 

= 10 4 , 

X* 

= 1, 

(HI) : 

AT* 

= 10 3 , 

X* 

= 10 2 


(87) 


( 88 ) 

(89) 

(90) 


all with /3 = 0.1 and which satisfy the small-scale constraints 
of Eqs. (l83l) and (l85l >. as well as the cosmological constraints. 
In particular, in each case we have 

(I) : Xs.s. - -10 6 , /3 2 er < 500, ^ - 1 ~ 10- 12 , (91) 

(II): Xs.s. ~ -10 4 , /3 2 er < 0.5, ^§-l~l(T 8 , (92) 

(III) : xs.s. - -10 6 , /3 2 er < 0.05, lO" 8 , (93) 

where Xs.s. is obtained from Eq. (l68l) and /3 per is the upper 
bound from Eq. (l70ll . Thus, /3 = 0.1 is consistent with the 
perihelion constraint for these three models. We can also note 
that the speed of scalar waves given by Eq.(l3Tl> is always very 
close to the speed of light (in the Solar System). Therefore, 
superluminality in the scalar sector is highly suppressed. 

We show these kinetic functions in Fig. |T] as a function of 
|x| - The models (I) and (II) correspond to a late transition 
from the K' ~ K* to K' ~ 1 regimes, with two possible val¬ 
ues for AT*. Therefore, during most of the cosmological evo¬ 
lution K' 1, which means that the scalar field is screened 
and departures from A-CDM are small. The model (III) cor¬ 
responds to an early transition, so that in a large range of red- 
shifts we have K' ~ 1 and deviations from A-CDM are of 
order /3 2 = ICC 2 . These behaviors are explicitly shown in the 
middle panel of Fig. U 

For simplicity we consider the simple rational function of 
Eq. (l82l) . which is even and leads to the same behavior in the 
two regimes x — > ±oo. However, there is a great freedom 
in the positive range x > 0- the only constraint being that 
K'(x) V' 1 for X > 1- For instance, we could add to Eq. (l86ll 
any function K+(x) that is negligible for (—x) A> 1 and sat¬ 
isfies K' + { 0) = 0, K'_ |_(x) + 3X* 1 at large x, such as 

K+(x) = exp(x 3 ). 

For comparison with previous works lfl3L(1714 1 911 . we also 
consider the purely cubic model 

(IV): K( X ) = ~l + X + KoX 3 - (94) 

This can also be seen as an effective model for moderate val¬ 
ues of x, probed by cosmology, while leaving the large nega¬ 
tive regime (—x) 1 unspecified, where the function K(x) 
needs to be modified as described above to satisfy Solar Sys¬ 
tem constraints. This model is also shown in Fig. Q] for the 
case Kq = 1. For most of the cosmologically relevant values 
of X> 0 < X < 10, the phenomenology of model (IV) lies 
between that of models (I) and (III). 


B. Theoretical consistency 

We have seen above that the models (I)-(III) satisfy the 
quantitative constraints associated with cosmological and So¬ 
lar System tests. In addition, they satisfy the generic theoret¬ 
ical consistency requirements. The four models (I)-(IV) have 
an even derivative K '{\;) > 0 and the functions W±(y) = 
yK'(±y 2 / 2) are monotonically increasing up to infinity over 
y > 0. As shown in Refs. E 0, this ensures that 
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these models are well behaved. More specifically, a static 
scalar field profile exists for any matter density profile (branch 
W- (y), see Eq. ( IT8l > ) and the background Klein-Gordon equa¬ 
tion can be integrated up to arbitrarily high redshifts where 
p — > oo. Moreover, there are also no ghosts nor small-scale 
instabilities. This corresponds to the constraints (|601)-(|6TT) and 
(l63li-(l64li. When the function K'(x) is even these two sets of 
constraints actually coincide. 

As noticed in Ea. (fTill . the speed of scalar waves around 
static backgrounds is greater than the speed of light. As seen 
from Eqs.(|9T1)-(|93]i. it is actually extremely close to c in the 
highly nonlinear regime, which applies to the Solar System. 
Superluminality is sometimes associated with possible viola¬ 
tions of causality, although such pathological behaviors can 
also be encountered within General Relativity (for exact so¬ 
lutions such as the Godel metric that do not describe realistic 
metrics), see the discussions in f27 :1 ] and |28|| . We argue in the 
appendix that space-time does not have closed time-like loops 
for realistic K-mouflage models, in particular for the models 
(I)-(III) of Egs. ( 1 8 8| >-(| 90 1). Therefore, the small superluminal¬ 
ity in the Solar System does not lead to causality problems. 


C. Cosmological evolution 


From the Klein-Gordon equation ([6|. one can show El 
that the time derivative of the cosmological background value 
(p of the scalar field scales with time as 




MpiK' Pt ’ 


and x 


p 2 H 2 


(95) 


At late times we have 

z^O: K'~ 1, x~/3 2 = 10- 2 , (96) 

^ ^m0 

where 0 \q and O ln o are the dark energy and matter cosmolog¬ 
ical parameters today. At early times we have, for the models 
(I), (II), and (III), 


z»l: K' ~ K,, x 


P 2 H 2 

kIh 2 ' 



(97) 



z 



z 



Thus, in the early matter era, the ratio p^/p goes to a small 
constant, which in our case is constrained by the Cassini 
bound Eq. (l66l > to be smaller than 10 -5 (this is because we 
assume the two regimes x ±oc have the same functional 
form). Thus, there remains a very small residual fraction of 
the scalar field energy density. For the model (IV) of Eq. ( l94l) . 
where K'(x ) g° es to infinity for x -A +oo, p v /p —> 0 at high 
redshift. In all cases, ptp/p ra d goes to zero at high redshift in 
the radiation-dominated era. 

At low redshifts, where \ falls below Xcan. = X*/v/7G, 
we have K’ ~ 1 and the model behaves like a canonically 
normalized scalar field coupled to matter in the presence of 
a cosmological constant. At higher redshifts, the nonlinear 
structure of the scalar field Lagrangian comes into play and 
K’ increases until it reaches K t . when x is of order x*. Then, 
the scalar field shows a K-mouflage-type screening, due to the 


z 


FIG. 2: Upper panel: variation of x{ z ) as a function of the redshift 
for the models of Fig. [T| All models have ft = 0.1 Notice the uni¬ 
versal behavior of x at small redshift given by d95l ) with K' ~ 1. 
Middle panel: variation of K' = d/\/dx(x) as a function of the 
redshift for the same models. Lower panel: relative deviation of the 
Hubble rate from the A-CDM prediction for the same models. 


large derivative prefactor K' that freezes the fluctuations of 
the scalar field. In terms of uniform background values, this 
screening still leads to the suppression of x and p v by factors 
l/K' 2 and 1/K'. 

The deviations from General Relativity and A-CDM pre¬ 
dictions are mostly sensitive to the ratio p 2 /K'. Therefore, 
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z 



z 




M [h A MJ 


FIG. 3: Upper panel: variation of ei{z) as a function of the red- 
shift for the models of Fig. Q] Lower panel: relative deviation of the 
growth factor from the A-CDM prediction for the same models. 


they should increase from models (II) to (I) to (III). Indeed, in 
models (I) and (II) the derivative K' quickly reaches the large 
value K t , because = 1, while in model (III) it remains 
of order unity over a large interval around \ = 0, because 
X* 1. The deviations from A-CDM are larger for model 
(I) than model (II) because of the smaller value of K, in the 
former. 

In Fig. [2] we plot the redshift evolution of \ and K', and the 
deviations of the Hubble rate from A-CDM. The figure shows 
that, as expected, the deviations from A-CDM increase from 
models (II) to (I) to (III). Also, in agreement with Fig. Q] the 
predictions of the cubic model (IV) lie in between the arctan 
models (I) and (III), as it gives an intermediate growth of K' 
as x increases at higher z. At 2 —> 0 all models give the same 
value for x, in agreement with Eq. ([95}, because x ~ /? 2 1 

and K' ~ K'( 0) = 1. The deviations from A-CDM for the 
Hubble rate peak around z ~ 1. All models are normalized 
to the same Hubble constant // () today and they also recover 
the same Einstein-de Sitter expansion at high redshift in the 
matter era, before dark energy becomes important. 

The upper panel of Fig.[3]shows that the four models satisfy 
the Lunar Ranging constraint on d In Q/At as the values of £2 
are below 0.01, in agreement with Eq. ( BTl) . This, however. 


FIG. 4: Upper panel: relative deviation of the power spectrum from 
A-CDM for the models of Fig. [T] at 2 = 0. Lower panel: relative 
deviation of the halo mass function for the same models, at z = 0. 


does not prevent the linear growth rate / = d In D+ /d In a 
today from being different in the different models, as shown 
in the lower panel. The deviations from A-CDM of the lin¬ 
ear fluctuations follow the same pattern as for the background 
cosmology. We also find that a maximal deviation of order 
2/3 2 is obtained for model (III) with a large value of x*, which 
behaves like a free scalar field coupled to matter in the recent 
past. The greatest deviations from A-CDM, at fixed /?, would 
be obtained for a model where K' ~ 1 over all relevant red- 
shifts. The results from such a model would be similar to 
those of model (III). 

In Fig. [4j we show the relative deviations from A-CDM of 
the nonlinear matter power spectrum and the halo mass func¬ 
tion. The method used for the computation of the power spec¬ 
trum was presented in ITbIi . It combines one-loop perturbation 
theory and a halo model and it has been tested against numer¬ 
ical simulations 0 of other modified gravity scenarios such 
as f(R), dilaton and symmetron models. The power spec¬ 
trum deviates on large scales by a constant which can reach 
a few percent. This reflects the boost in the linear growth of 
structure depicted in the lower panel of Fig. [3] In the non¬ 
linear regime, mode coupling helps to boost the deviations 
from standard A-CDM even further. In the figure, the rela- 
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tive difference starts to decrease after k ~ 1 h/Mpc, which 
is a consequence of having used the same halo concentration 
parameters in the halo model formalism 0. In reality, the 
K-mouflage field should also affect the concentration of dark 
matter haloes, which should translate in modifications of the 
small-scale clustering power as well. The enhanced gravita¬ 
tional strength should lead to a deepening of the gravitational 
potentials of the haloes, making them therefore more concen¬ 
trated (see e.g. |[3(ll3lll for an example of this in the context of 
Galileon and Nonlocal models of gravity, respectively). This 
would increase the relative difference of Pj. at high -k values. 
On the other hand, the enhanced forces may also make parti¬ 
cles inside haloes move faster. This can cause the halo to ex¬ 
pand slightly, and therefore, to become less concentrated (see 
e.g. 122] for an example of this in coupled quintessence mod¬ 
els, which are similar to unscreened K-mouflage models). A 
detailed investigation of this effect involves running dedicated 
A’-body simulations which is left for future work. 

For the halo mass function, we take into account the_ef- 
fect of K-mouflage through the spherical collapse model |18|. 
The amplification of gravity by the fifth force implies that a 
smaller initial density fluctuation Sli at a high redshift z\ is 
needed to produce a collapsed halo of a given mass M at low 
z, as compared with the A-CDM reference. Moreover, such 
objects are less rare because of the faster growth of structures 
(e.g., of the matter power spectrum). Both effects amplify 
the large-mass tail of the halo mass function, as seen in the 
lower panel of Fig. 0] [the mass function is written in terms 
of the usual scaling variable v = with the char¬ 

acteristic Gaussian tail n{M) ~ e~ v / 2 at large mass, and 
the fifth force implies both a smaller critical threshold 8l and 
a greater rms density fluctuation cr(M)]. The mass function 
is slightly decreased at low masses because of mass conser¬ 
vation, as the integral over the halo mass function is normal¬ 
ized so as to give back the mean density of the Universe. The 
lower panel of Fig. 0] shows that K-mouflage models can ex¬ 
hibit a 5 — 20% boost in the number density of cluster mass 
haloes, 10 14 M Q /h < M < 10 15 M Q /h. This may offer a 
clear enough signal to be observable with future surveys. 


VIII. CONCLUSION 

We have determined the conditions for K-mouflage models 
of gravity to satisfy the stringent Solar System tests of gravity, 
whilst remaining sufficiently different from standard A-CDM, 
and hence, cosmologically interesting. In particular, we have 
used the results from Cassini on the amplitude of fifth forces in 
the Solar System, and bounds from Lunar Laser Ranging ex¬ 
periments on the anomalous perihelion and the rate of change 
of the gravitational strength in the Earth-Moon system. 

We showed that the conformal coupling strength of the 
K-mouflage field to matter is constrained as /3 < 0.1 
(cf. Eq. d58t ). by the condition that the time variation of New¬ 
ton’s gravitational strength in the Jordan frame is compatible 
with the Lunar Ranging bounds. This constraint is indepen¬ 
dent of the detailed functional form of the K-mouflage La- 
grangian density term K(x) [cf. Eq. ©], and is also tighter 


than the bounds coming from BBN (cf. Eq. (l55l >). By focus¬ 
ing on static configurations, which correspond to the branch 
X < 0 of the function K (x), we have seen that the perihelion 
bound is more stringent than the Cassini result in constraining 
the functional form of K(x)- hi particular, for K-mouflage 
models to remain compatible with these tests, any nonlinear 
terms in K(x < 0) should be highly suppressed (cf. Eq. (1701 )). 
For instance, models for which K (x) ex x 3 will fail to meet 
the Solar System requirements. On the other hand, the shape 
of the coupling function A(ip) is not strongly constrained. 

We have presented several explicit models that satisfy these 
Solar System and cosmological bounds as well as generic the¬ 
oretical consistency conditions. In particular, these models 
have well-defined solutions up to arbitrarily high redshift, in 
the cosmological context, and for any matter density profile, 
in the small-scale context. There are no ghosts nor small-scale 
instabilities. Even though scalar waves can propagate at a 
speed that is slightly greater than the speed of light around 
small-scale static backgrounds (with a relative difference of 
only 10 -8 or less in the Solar System), there are no closed 
time-like loops nor causality problems, from the Solar Sys¬ 
tem to cosmological scales. 

One of our main results is that, despite the constraints on /? 
and K (x) from the Solar System tests, one is still able to find 
a family of K(x) functions (cf.Eq. (f82V ) that has interesting 
and potentially testable cosmological predictions. We have in¬ 
vestigated the main cosmological aspects of the models char¬ 
acterized by Eq. (l82l >. For a set of illustrative cosmological 
parameters, we have seen that the linear growth of large-scale 
structures can be boosted by a few percent by the present day 
(cf. Figs.0]and0]i. Our results from semi-analytical models of 
structure formation also show that this difference gets ampli¬ 
fied on smaller scales, where the evolution of the matter den¬ 
sity field becomes nonlinear. Moreover, the expected number 
density of cluster mass haloes shows also a 5 — 20% enhance¬ 
ment, relative to A-CDM (cf. Fig. 0}. Another interesting 
aspect of these models is that their expansion history can be 
slightly different from the A-CDM scenario, with deviations 
at the percent level or slightly below that may be constrained 
by observations (cf. Fig. 01). This is different from the cases of 
DGP and/or f(R) models of gravity, where the expansion can 
follow the A-CDM scenario up to very high accuracy or even 
exactly. This means that the parameter space of these models 
can be constrained by the position of the acoustic peaks of the 
CMB temperature power spectrum, as investigated already in 
0, but for models that fail Solar System tests. 

To summarize, the models we have built up are predic¬ 
tive and distinguishable from other alternatives to A-CDM, in 
the perspective of future experiments such as Euclid 0 and 
LSST [211]. In the future, 21-cm intensity mapping both dur¬ 
ing [33] and after the completion of the reionization 1040361] 
will open new windows to test modified gravity and will help 
in discriminating between models and in constraining further 
the shape of the K-mouflage function K(x). Future 21-cm 
surveys such as with the Square Kilometre Array will probe 
the Universe’s expansion up to higher redshifts and the mat¬ 
ter power spectrum down to smaller scales, especially in the 
range 2 < z < 8 of interest for the K-mouflage model. 
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In future work, we believe it would be of interest to per¬ 
form more focused studies of the cosmological constraints in 
these models, by following, for instance, the line of work of 
0. It would also be interesting to study more accurately 
the predictions for nonlinear structure formation by running 
A ? -body simulations. Such studies should provide a clearer 
picture of how these types of modifications to gravity can im¬ 
pact on several cosmological observables, which should help 
in the interpretation of the results from future observational 
missions. 
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Appendix A: Superluminality and Causality 


In this appendix, we elaborate on the superluminality of 
scalar perturbations and its link with causality. We have 
seen in the main text that K-mouflage models passing the 
Solar System tests are such that scalar perturbations around 
a static background propagate with a speed greater than the 
speed of light. This may cause instabilities and in particu¬ 
lar a loss of causality with signals being transmitted to the 
past, i.e. the existence of time-like closed curves. This issue 
was tackled for K-essence models in 0. We follow a sim¬ 
ilar method here. Let us first expand the K-mouflage action 
to second order in n, where ip = (p + n. Here ip(x, t) is a 
background configuration that may depend on scale and time 
[in the small-scale static regime f ix) only depends on posi¬ 
tion, whereas in the large-scale cosmological regime (p{t) is 
the homogeneous cosmological background]. We also denote 

X = -g^d^(pd v (pl2M A , K' = K’(x) and K" = K"(x). 

The second order part of the action reads 


Vo = 


d 4 xy/—g 


-K'd^ndfj, 7r + 7r) 2 


(Al) 


It is convenient to define the disformal metric 


The disformal metric G t , v is the metric felt by the scalar per¬ 
turbation and we can check from Ea. (lA4l) that it is Lorentzian. 
Therefore, initial-value problems for it are well posed on any 
smooth space-like Cauchy surface E for the metric G t , u , and 
the solution is unique and propagates causally (see Sec. 10 in 
I2l). In the static case, we retrieve that the Klein-Gordon 
equation is hyperbolic with a propagation speed c s given by 
Eg. (OH (for small wavelengths), whereas around the cosmo¬ 
logical background we recover the propagation speed c s given 
by Ea. (l59l > (for high frequencies). 

In general, the propagation of n occurs in the disformal 
metric G IJ V . Space-time equipped with the metric G IUJ is sta¬ 
bly causal, i.e. there are no time-like closed loops (including 
for infinitesimal deviations from the metric provided 

there exists a globally defined function / on all space-time 
which is time-like, i.e. G^d^fd^f < 0 10 . Following 
10 , we look for a “global time” / that applies to both geome¬ 
tries g llu and G II V and thus guarantees the absence of closed 
causal loops. A simple choice is to choose the cosmic time 
t, which clearly satisfies the required property for the metric 
g Considering the Newtonian gauge, which describes all 
systems that we study in this paper, from the cosmological 
background and perturbative regime down to the Solar Sys¬ 
tem, 


ds 2 = g^dx^dx" = — (1 + 2^<^)dt 2 + a 2 (i)(l — 2H/N)dx 2 , 

(A6) 

where T'n is the Newtonian potential, we have 

g^d^tdj = 1 < 0 for H7 N > -1/2. (A7) 

1 + 2Wn 


Since we focus on systems with |4 /n| -C 1 (e.g., H/n ~ 10~ 6 
in the Solar System), we have g^d^td^t < 0. On the other 
hand, we obtain 


G^drfdvt = - 


K'{1 + 2’J'n) + K" (doip) 2 / M a 
7(1 + 2T' n ) 2 


(A8) 


G' J ' V = 7- 1 (KY V - < A2 > 

with 

7 = (K'f/ 2 (K' + 2xK") 1/2 > 0. (A3) 


whence 

G^dJdG < 0 for C = K’ + K" ^f > 0, (A9) 

Ad 4 

where we used the approximation 1 + 24 /n ~ 1. Around the 
cosmological background, where x = (dip/dt) 2 /2A4 4 , we 
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obtain C = K' + 2 \K", whence C > 0. Around a static 
background, we obtain C = K' whence C > 0. 

For more general backgrounds, we can see from Ea. (lA9t 
that C > 0 as soon as K" > 0, which for the models (I)-(IV) 
of Eas. (l88l >-(|90} and ( l94l ) corresponds to \ > 0. On the semi¬ 
axis x < 0, we have seen that C ~ K' > 0 in the static limit, 
\dCp/dt\ < |V<)?|. Therefore the remaining case corresponds 
to \dCp/dt\ ~ |Vj3|, where time and spatial derivatives are of 
the same order, and with x < 0, that is, \d(p/dt\ < |V£|. 
Then, we have 

X < 0, K" < 0: C>K'-\K"x\, (A10) 

as we assume \ ~ —(d(p/dt) 2 /2M A ~ — (V<j?) 2 /2A4 4 . In 
the linear unscreened regime, |x| <C 1, this gives C > K' ~ 1, 
whence C > 0. In the highly nonlinear screening regime, 
|x| 2> 1 [more precisely x < — X* as i n <l65l) 1, we have 
seen in Eqs. (f75l i and ([79} that \xK"\ -C K' , because of the 
perihelion constraint. Therefore, on these nonlinear scales 
we obtain C ~ K' > 0, whether we are in the static limit. 


\d(p/dt\ <C |Vj5|, or not, \d<p/dt\ ~ |V^|, and we have al¬ 
ready seen that C > 0 when \d(p/dt\ > |V<^| because it im¬ 
plies x > 0 and K" > 0 in our models. Note that the regime 
\d(p/dt\ ~ |V(^| is unlikely to occur in practice in small- 
scale systems, because the quasistatic approximation applies 
very well, even for relatively fast matter density evolutions 
with matter flow velocities of order v ~ c/10 [13]. There¬ 
fore, we usually have \dCp/dt\ -C |V(^|, which directly gives 
C ~ K' > 0. 

Thus, we conclude that g IJU d^td v t < 0 and G^ u d^td v t < 
0 and there are no closed causal loops around usual astrophys- 
ical and cosmological backgrounds with > —1/2. This 
analysis fails close to neutron stars or black holes, where 'l'\- 
becomes large, but this is not related to the K-mouflage model 
as it already appears in the metric g in Eq.( lA7l ). Then, one 
must look for another global time coordinate, or over a large 
volume around the compact object, but we leave this analysis 
of more extreme astrophysical situations to future work. 
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